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Abstract 

We justify the Hamilton least action principle for the Maxwell-Lorentz equations with 
Abraham's rotating extended electron. The main novelty in the proof is application of the 
variational Poincare equations on the Lie group SO (3). The variational equations allow to 
derive the corresponding conservation laws from general Nother theory of invariants. 
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1 Introduction 



We justify the Hamilton least action principle for the system of Maxwell-Lorentz equations with 
a rotating charged particle. Our main contribution is the variational derivation of the Lorentz 
torque equation, see equation (jl.3p below. 

First recall the case of a finite system of material points (g^mj). The angular momentum is 
defined by 

M := J^Qi A Pi : = Am i<li- (1-1) 

By the second and the third Newton laws this implies 

M = J2<H Aft = J> A Ft = J2<H AF i Xt = T, (1.2) 

where T is called the external force torque. Our aim is to derive the similar torque equation for a 
charged rigid body in the Maxwell field: 

Iu = e J(x-q)A[E + E ext + (q + u A (x - q)) A (B + B ext )] p(x - q)dx, (1.3) 

where I is the moment of inertia, u is the angular velocity, and p(x) is a charge distribution, and 
the right hand side is the torque of the Lorentz force. 

Formally the rigid body can be considered as an infinite system of material points. However 
equation (11.31) cannot be obtained directly from (11.21) . since we cannot correctly take into account 
all the forces of mutual interaction between the different pieces of the rigid body. That is why 
we look for a different approach to the derivation of (ll.3p . We show that (II. 3p follows from the 
Hamilton variational least action principle with the standard interaction term —A p + A- j in the 
Lagrangian density. 

For the free rigid body (E, E ext , B, B ext = 0) equation (jl.3p reduces to the Euler's equations 
which have been obtained from the variational principle first by Poincare [10], and in [2] for an 
external force field with an axial symmetry. 

Note that in our case the fields E and B are generated by the motion of the charged body, 
and E ext , B ext are external fields. 

Let us comment on previous works. Equation ( II .3p is well recognized since the Abraham's 
works [3j H]. In [3j Section 11] Abraham computed the Lagrangian as integral of —A p + A - j for 
standing rotating spherically symmetric electron subject to external fields obeying very special 
symmetry conditions. In this case the Lagrangian depends only on one variable u), the angular 
velocity. However, derivation of the torque equation (II. 3p from the variational it Hamilton least 
action principle remained an open question. 

Nodvik applied the variational Hamilton principle to to the rotating charge in the Euler angles 
He deduced the dynamical equation [HI (5.46)] for co(t) and have established the corresponding 
conservation laws. However, the Nodvik equation looks differently from (ll.3p . Appel and Kiessling 
[5] suggested a transformation of [9, (5.46)] to (11. 3p . but they did not give the details, see [5j A. 1.4]. 
The direct derivation of the conservation laws from ( II .3p is presented by Kiessling in [8]. 

We propose an invariant derivation of (II. 3p from the Hamilton least action principle relying on 
Poincare equations [10112] on the Lie group SO (3). We also deduce the corresponding conservation 
laws by general Lagrangian formalism using the Nother theory of invariants. 

The plan of our article is as follows. In Section 2 we state the Maxwell-Lorentz equations for the 
rotating charge, and introduce the corresponding Lagrangian functional. Further we deduce the 
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equations from the Hamilton variational principle relying on the Poincare theory [10] presented in 
the Appendix A. In Section 3 we prove the conservation laws. Finally, in Appendix B we identify 
the commutators of invariant vector fields. 

2 Maxwell- Lorentz equations 

The Maxwell field consists of the electric field E(x,t) and the magnetic field B(x,t) generated 
by a motion of a rotating charge. The external fields E ext and B ext are also generated by the 
corresponding external charges and currents. Let the rotating charge be centered at the position 
q with the velocity q. For simplicity we assume that the mass distribution, m p(x), and the charge 
distribution, e p(x), are proportional to each other. Here m is the total mass, e is the total charge, 
and we use a system of units such that m = 1, e = 1. The coupling function p(x) is a sufficiently 
smooth radially symmetric function of fast decay as \x\ — > oo, 

P(X) = Pr (\x\). (C) 

2.1 Angular velocity 

Let us denote by u(t) E IR 3 the angular velocity "in space" (in the terminology of [2]) of the 
charge. Namely, let us fix a "center" point O of the rigid body. Then the trajectory of each fixed 
point of the body is described by 

x(t) = q{t) + R(t)(x(0)-q(0)), 

where q(t) is the position of O at the time t, and R(t) E SO(3) . Respectively, the velocity reads 

x(t) = q(t) + R(t)(x(0)-q(0)) = q(t) + R(t)R~ l (t)(x(t) - q(t)) = q(t)+u(t) A (x(t)-q(t)), (2.1) 

where u(t) E IR 3 correponds to the skew-symmetric matrix i?(t)i? _1 (t) by the rule 

/ -w 3 (t) u 2 (t) \ 
R{t)R-\t) = Ju(t) := w 3 (t) -ui{t) . (2.2) 

\ -u 2 {t) / 

We assume that x and q refer to a certain Euclidean coordinate system in IR 3 , and the vector prod- 
uct A is defined in this system by standard formulas. The identification (12.21) of a skew-symmetric 
matrix and the corresponding angular velocity vector is true in any Euclidean coordinate system 
of the same orientation as the initial one. 

2.2 Dynamical equations 

Then the system of Maxwell-Lorentz equations with spin reads, see [TT] 

E = V AB-(q + uA(x-q))p{x-q), B = —V A E, (2.3) 

V ■ E(x,t)= p(x-q(t)), V ■ B(x, t) = 0, (2.4) 

q = J [E + E ext + (q + u A (x - q)) A (B + B ext )]p(x - q)dx, (2.5) 
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Iu= J {x - q) A [E + E ext + (q + uA(x-q))A(B + B ext )]p{x - q)dx, (2.6) 
where / is the moment of inertia defined by 

1 = x 2 p(x)dx. (2.7) 

Here the equations (12.31) are Maxwell equations with the corresponding charge density and current, 
equations (12. 4p are constraints. The back reaction of the field onto the particle is given through the 
Lorentz force equation (12.51) . and the Lorentz torque equation (12.61) deals with rotational degrees 
of freedom. 

2.3 Lagrangian functional and variational principle 

Our main goal is to deduce equations (I2.3p - (I2.6I) from the Hamilton least action principle. First 
let us introduce electromagnetic potentials A = (Aq,A), A ext = (Aq**, A ext ): 

B = V A A, E = -VA, - A. (2.8) 

B ext = V A A ext , E ext = -VA ext - A ext . (2.9) 
Next we define the Lagrangian 

L(A, q, R, A, q, R) = l -J (E 2 - B 2 ) dx + l -q 2 + hu? - J [A Q + A e xt ]p(x - q)dx + 

j \q + u A(x-q))-[A + A ext }p{x - q)dx, (2.10) 

where E, B are expressed in terms of A, A by ( \2.8L and uj = J^RR~ l by ( 12.2ft . 
The last two integrals represent the interaction term 

j[(A + A e xt }p-j[A + A ext }dx 

in view of ( 12.11) . The corresponding action functional has the form 



S = S(A, q, R) := L{A{t),q{t),R{t),A{t),q{t),R{t))dt (2.11) 



Then the Hamilton least action principle reads 

5S(A,q,R) = 0, (2.12) 

where the variation is taken over A(t),q(t), R(t) with the boundary conditions 

(5A,5q,5R)\ t=tl = (5A,8q,5R)\ t=t2 = 0. (2.13) 

We assume that all the involved functions and fields are sufficiently smooth and have (with all 
the necessary derivatives) a sufficient decay as |x| — > oo so that partial integrations below could 
be possible. 

Our main result is the following theorem. 
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Theorem 2.1 The Maxwell- Lorentz system with spin A2.3\) to A2.6\) is equivalent to the least 
action principle A2.12\) - f2.13\) . 

We will analyze the variations in A, q, R separately, namely, we prove that 

u = ° <•>■ f = ° < 6) - i = ° w (2 - 14 > 

is equivalent to f l23|) -f l23|) . 

2.4 Equations for fields and particle trajectory 

Equations for fields It is well-known [6] that ( 12.14ft (a) is equivalent to 

d SL . , 

^77 = L * 2 ' 15 
at 8 A 

and that the last Euler-Lagrange equations are equivalent to the Maxwell equations (I2.3P with the 
constraints (12. 4p . 

Remark Note that the terms with Aq x1 and A ext in (12.101) are additive and remain additive 
while one makes variations in q, and R. Then, for simplicity of exposition, we put in all of the 
computations below A^ xt = and A ext = 0. 

Equations for particle trajectory Similarly, (12. 14ft (b) is equivalent to 

jL, = L,. (2.16) 

It remains to check that the last Euler-Lagrange equations are equivalent to the Lorentz force 
equation ( 12.51) . Let us change variables in the last two integrals of ( 12.101) and obtain 

Aq(x + q, t)p(x) + J (q + u A x) ■ A(x + q, t)p(x)dx. 

We have 

Lg = q + J A(x + q,t)p(x)dx, L g = — J \7A (x + q,t)p(x)dx + J(q + uAx) ■ V A(x + q,t)p(x)dx. 

Here, for a vector field b(x) we denote b ■ VA = J^bjVAj. Then the Euler-Lagrange equation 
( 12~T6|) reads 

q+ J (^A(x + q,t) + (q- V)A(x + q, t)j pdx = - J VA (x+g, t)p(x)dx+ J (q+uAx)-V A(x+q, t)pdx. 
Substituting A(x + q,t) = —E(x + q,t) — VA (x + q, t) we obtain 

g = y E(x+q,t)p(x)dx+ J (q-VA(x + q,t)-(q-V)A(x + q,t))pdx+ J(uAx)- VA(x+q, t)pdx. 

' (2.17) 
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To make the notations shorter let us omit the dependence on t in the further computations. By 
the identity q ■ VA(x + q) — (q ■ V)A(x + q) = q A V A A(x + q) we obtain 

j \q-VA(x + q) - q-VA(x + q))pdx = J q A V A A(x + q)pdx. (2.18) 

It remains to check that 

(co A x) ■ VA(i + q)pdx = IwAiAVA A(x + q)pdx. (2.19) 



Let us check for the first component, for the rest the computation is similar. The first component 
of the LHS of fl2~T9l) equals 

J [(u 2 x 3 - u 3 x 2 )d 1 A 1 (x + q) + (u 3 x 1 - u 1 x 3 )d 1 A 2 (x + q) + {ijO\x 2 - u 2 x 1 )d 1 A 3 (x + q)) pdx. 

The first component of the RHS of (I2.19P equals 

[{uj 3 xi - uix 3 )(d 1 A 2 (x + q) - d 2 A 1 (x + q)) - (uoix 2 - u 2 x 1 )(d 3 A 1 (x + q) — diA 3 (x + q))] pdx. 

For the difference of the LHS and the RHS we apply partial integration, and obtain 

[(u 2 x 3 - u 3 x 2 )diAi(x + q) + (w 3 £i - uix 3 )d 2 Ai(x + q) + {uix 2 - u 2 xi)d 3 A 1 (x + q)] pdx = 

-J Ai(x + q) [ui{x 2 d 3 - x 3 d 2 ) + u 2 {x 3 di - xid 3 ) + Uz{xxck - x 2 d{)} pdx = 

- J ' Ai(x + q)(u ■ V e )pdx, 

where Vg = (V^, Ve 2 , Ve 3 ) and V# . is the differentiation in the angular coordinate 9j around the 
coordinate axis Xj. Since p is radially symmetric, one has VqP = and we come to (I2.19p . From 
(I2T81 . fl2TT8|) . and fl2~T9l) we conclude that the equation (12T7I) reads fl2~5l) . 

2.5 The torque equation 

Finally, it remains to check that (12 .14ft (c) is equivalent to (12.61) . 

To make the corresponding variation in R, let us express u(t) = J _1 i?(t)i?~ 1 (t) in the right- 
invariant vector fields on SO(3). Namely, consider an orthonormal basis {e&} with the right 
orientation in IR 3 . Then 

eiAe 2 = e 3 , e 3 A e x = e 2 , e 2 Ae 3 = e 1 . (2.20) 

Let us express the angular velocity in {e^}: u(t) = J2uk(t)ek- The space so(3) of skew-symmetric 
3x3 matrices with the matrix commutator is isomorphic to IR 3 with vector product by the 
isomorphism J of (12.21) : 

—LU 3 UJ 2 

u 3 -ui I = J(ui,u 2 , co 3 ). (2.21) 
— u 2 uji 
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In detail, if A, B e so(3), a, b e H 3 , and A = Ja, B = Jb by the isomorphism (I2.2ip . then 

AB — BA = J{a A b). (2.22) 

Further, R{t)R~ l {t) e TeSO(3) is the tangent vector R{t) of 50(3) at the point R{t) translated 
to the unit E of S0(3) by the right translation i? _1 (t). By the linear isomorphism (12.211) . 

R^R-^t) = J2Mt)e k , e k := J~ x e k . (2.23) 

Then 

R(t) = R(t)R-\t)R(t) = u k (t)v k (R(t)), v k (R) := ~e k R. (2.24) 

As the result, R(t) has the same coordinates w.r.t. the vector fields v k at the point R(t) as 
uj(t) in the basis {e^}. The fields v k {R) are right translations of e k and hence are right-invariant. 
The next lemma is proved in Appendix B. 

Lemma 2.2 For the constructed above vector fields v k on SO (3) the following commutation rela- 
tions hold: 

[v 1 ,v 2 ] = -v 3 , [v 3 ,vi] = -v 2 , [V2,v 3 ] = -vi. (2.25) 

Poincare equations Now we are going to deduce (12.61) from the Poincare equations. Namely, 
as shown by Poincare [TOl [2] , the equation (12.141) (c) is equivalent to the Poincare equations 

ij J 

where L is L with u(t) expressed in the coordinates (u)i(t), a; 2 (t), 0J 3 (t)). For convenience of the 
reader the derivation of Poincare equations is presented in Appendix A. 

Note that the Lagrangian L does not depend explicitly on R, and hence v k {L) = 0, k = 1, 2, 3 
by (I A .40 . Hence, the corresponding Poincare equation reads 

d dL x - j dL 

jtd^ = X? aUi d^' 



Equivalence to the torque equation It remains to check that the equations (I2.26P are equiva- 
lent to the Lorentz torque equation (12. 6p . It suffices to check the equivalence for the first component 
with k = 1, since for the rest k the computation is similar. 



First, 

dL 



Therefore, 



Iuj x + / (x 2 A 3 (x + q) - x 3 A 2 (x + q))p(x)dx. 



d dL 
dt du\ 



[x 2 (A 3 (x + q) + q- VA 3 (x + q)) - x 3 (A 2 (x + q) + q ■ VA 2 (x + q)))pdx. 
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Second, by (12 . 25[) and (1A.5|) one has c 3l = —1, = 1, and all the rest c^=0. Hence, 



£4 



dL 



dL dL 



0J2- 



duj du 3 



u 3 



duo 



Iu 3 + j ' (x 1 A 2 (x + q) - x 2 Ai(x + q))pdx — oj 3 Iu 3 + J (x 3 Ai(x + q) - XiA 3 (x + q))pdx 

J [-{u 2 x 2 + uj 3 x 3 )Ai(x + q) + xi(u 2 A 2 (x + q) + u?,A 3 (x + g))] pdx. 
Finally, we come to the equation 

Iui = / (x 3 A 2 (x + q) - x 2 A 3 (x + q))pdx + / [x 3 (q ■ V)A 2 (x + q) — x 2 (q ■ V)A 3 (x + q)]pdx+ 



J [xi(u 2 A 2 {x + q) + u 3 A 3 (x + q)) - (u 2 x 2 + u 3 x 3 )Ai(x + q)) pdx. (2.27) 
Now let us proceed to the first component of the equation (12.61) . Using the identity 

x A [(u A x) A B\ = (u A x)(x ■ B) 

we obtain 

Iu 1 = J (x A E(x + q))ipdx + J (x A (q A B(x + q)))xpdx + J {u Ax) x (x ■ B(x + q))pdx, (2.28) 
We insert 5 = V A A and obtain that the first integral equals 



x 3 A 2 (x + q) - x 2 A 3 (x + q) + {x 3 d 2 - x 2 d 3 )<p 



pdx 



j ' (x 3 A 2 (x + q) — x 2 A 3 (x + q))pdx - J (p(x 3 d 2 - x 2 d 3 )pdx = 

j ' (x 3 A 2 (x + q) - x 2 A 3 (x + q))pdx, 
since (x 3 d 2 — x 2 d 3 )p = do 1 p = 0. Further, the second integral of (12.281) reads 

J [x 2 {qi{d 3 A x {x + q)- diA 3 {x + q)) - q 2 (d 2 A 3 (x + q) - d 3 A 2 (x + q))- 

x 3 (q 3 (d 2 A 3 (x + q)- d 3 A 2 (x + q)) - q^A^x + q) - d 2 A x (x + q))}pdx. 
Finally, the third integral of (12.281) gives 

j [xi(u 2 A 2 (x + q)+ cu 3 A 3 (x + q)) - (u 2 x 2 + w 3 x 3 )Ai(x + q)]pdx+ 



(2.29) 



(2.30) 



(u 2 x 3 - u 3 x 2 )(Ai(x + q)(x 3 d 2 - x 2 d 3 ) + A 2 (x + q)(x x d 3 - x 3 d x ) + A 3 (x + q){x 2 d x -x x d 2 ))pdx 



J [xi{u 2 A 2 (x + q) + u 3 A 3 (x + q)) - (u 2 x 2 + u 3 x 3 )Ai(x + q)]pdx, (2.31) 

since 

{A x {x + q)(x 3 d 2 - x 2 d 3 ) + A 2 (x + q)(x 1 d 3 - x 3 d x ) + A 3 (x + q)(x 2 d 1 - xid 2 ))p = A(x + q) ■ W e p = 0. 
By (12351) . (CTD , ([MID, the difference between the RHS of QX27p and the RHS of fl2T2gj) equals 

J [x 2 (gi5 3 v4i(x + g) + g 2 <9 3 v4 2 (x + q)) - x 3 {q 3 d 2 A 3 {x + q) + gi<9 2 ,4i(a; + g))] pefe+ 

/ [x 2 g 3 9 3 A 3 (x + g) - x 3 q 2 d 2 A 2 (x + g)] prfx = 

j q ■ A(x + q){x 3 d 2 - x 2 d 3 )pdx = 
and we obtain that the equation (12.271) reads (I2.28p . The theorem is proved. □ 

3 Conservation laws 

We have derived the system fl2.3l) - fl2.6p by the least action principle with the Lagrangian fl2.10p . 
When the external fields posess a symmetry, the correponding conservation laws can be also 
derived from the Lagrangian formalism. Let us recall that the Lagrangian (I2.10p reads 

L(A, q, R, Aq,R) = \J (E 2 - B 2 ) dx + + hu 2 - 

J {A + A e xt )p(x - q)dx + J{q + u A (x - q)) ■ (A + A ext )p{x - q)dx, (3.1) 

where u = J^RR' 1 , A = (A , A), A = (A , A), and 

E = -VA Q - A, B = VAA. (3.2) 

As above, we denote 

L(A, q, R, A, q, u) = L(A, q, R, A, q, R) (3.3) 
where u = (ui,co 2 ,u 3 ) is defined by (12.241) i.e. Uk are coordinates of R in the basis Vi(R),v 2 (R), v 3 (R). 

3.1 Energy conservation 

Let us denote 

X := (A,q,R), V := X= (A,q,R). (3.4) 

Let Ay xt end A ext do not depend on time. Then the Lagrangian (13.11) does not depend on t, and 
the energy 

E(X, V) := L v - V — L (3.5) 
is conserved, pQ. By (13. 4 j) and since L does not depend on A , we have 

L v ■ V = L A ■ A + Lg ■ q + L k ■ R = L A ■ A + Lg ■ q + L w ■ u. (3.6) 
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Proposition 3.1 The energy reads 

E= l -J {\E\ 2 + \B\ 2 )dx + + hu 2 + J A e xt p{x - q)dx. (3.7) 
Proof By CI3 .31) and (13.11) . one has 

L A - A = - J E ■ Adx, Lg-q = q 2 + / q-(A + A ext )p(x - q)dx, 



and 



L u ■ u = Iu 2 + J co A(x-q) ■ (A+ A ext )p(x - q)dx. 



Then 

E = L^-A + Lq-q + L UJ -u-L 

= X -q 2 + hu 2 + ~ y (|5| 2 - |£| 2 )dx + J(-E-A + A -{V-E) + A e xt p{x - g))dx. (3.8) 
The last integral equals 

- J{E-A + VA • £ + A ext p{x - q))dx = - J (E(A + VA ) + A ext p{x - q))dx 

= J E 2 dx + J A e xt p(x - q)dx 
and hence (I3.8P reads (13. 7p . □ 

3.2 Momentum conservation 

Let us consider the spatial translations of the lagrangian coordinate X = (A, q, R): 

(A(x), q, R) \-> (A(x - h),q + h,R). 

If the external field A ext = (AQ Xt ,A ext ) does not depend on Xj with some j, then the Lagrangian 
(13.11) is invariant w.r.t to the one-parametric group of spatial translations 

gi(A(x),q, R) = (A(x — sej),q + sej, R), (3.9) 

where ej G IR 3 is the corresponding basis vector. By the Nother theorem pQ the expression 

Pi = PjiX, V) := L v ■ ^P-\s=o (3.10) 

is conserved. 

Definition 3.2 Vector P = (Pi, P2, P3) is called momentum of the state (X, V). 
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Proposition 3.3 The momentum reads (cf. f^j) 

P = q + / EABdx+ / A ext p(x-q)dx. (3.11) 



Proof For concreteness let us compute Pi. Formula ( I3.9P implies 



\s=0 



-(ei-VA(ar), e x ,0). 



ds 

Since L does not depend on A , and the map g\ leaves u unchanged, 

dgXX), 

Pi = L v - -^^| s =o = -L A ■ (ei ■ V)A + ■ ei = 
- y (VA + A) ■ (ei • V)Adx + gei + ^ e x • - g)dx + J A\ xt p(x - 

Note that 

qi+ Aip(x - q)dx - \ (VA + A) ■ Q X A dx = 



qi+ J A lP (x-q)dx-J (0 1 Ao0 1 A 1 +0 2 A 1 A 2 +0 3 A 1 A 3 +A 1 1 A 1 +A 2 1 A 2 +A 3 1 A 3 )dx. (3.12) 
On the other hand, consider the RHS of (13.111) and obtain 

Qi + J (E A B^dx = Ql + J [(-d 2 A - A 2 )(d 1 A 2 - d 2 A x ) + (0 3 A + A^idsAt - d x A 3 )]dx = 



qi + J (-d 2 A d l A 2 -A 2 d 1 A2 + d 2 A d 2 A 1 +A 2 d 1 A2+d 3 A d 3 A 1 +A- i d 3 A 1 -d 3 A d 1 A- i -A 3 d 1 A 3 )dx. 

(3.13) 

The difference between ( 13.12ft and ( 13.131) equals 

J A lP {x-q)dx- J (d^od^ + 0^08^ + 0^083 A, + A 1 d 1 A 1 + A 2 d 2 A 1 + A 3 d 3 A 1 )dx. (3.14) 
Finally, the partial integration implies 

A 1 p(x-q)dx= / A 1 (V-E)dx= / A x V • (-VA - A)dx = I Ai(-AA - VA)dx = 



J {d\Aod\Ai + 2 A 2 A x + 3 Ao0 3 A x + A.O^ + A 2 2 A x + A 3 3 A x )dx 

Hence, the difference (13. 14ft equals zero and the first components of the LHS and RHS of (13.111) 
are equal. □ 
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3.3 Angular momentum conservation 

Let the external potential A ext be axially symmetric, 

A ex \U k x) = U k A ext (x), (3.15) 
where U k is any rotation around the axis x k . 

Lemma 3.4 Let Ii3.15\) hold. Then the Lagrangian Ii3. 1\) is invariant w.r.t. the axial rotations 

A {x) ^ MU^x), A(x)^U k A(U l ; 1 x) 1 A(x) m- U k A(U^x); (3.16) 

q i— >■ U k q, q i— >■ f/ fc g, (3.17) 
R^U k R, R^U k R, (3.18) 

Proof By (13.21) the transforms (13. 16[) of the potentials induce the following transforms of the 
fields: 

E(x) ^ U k E{U^ l x), B(x) ^ U k B{U k l x). (3.19) 

Further, we have, in operator notations, Juo = uA, where wA is the operator of the vector product 
by u in IR 3 . Then it is easy to check that J{U k uj) = U k J'(co)U k ~ l . Thus, for u = J^RR^ 1 we 
obtain Ju = RR' 1 and hence J{U k u) = U k {RR- l )U k l = {U k R){U k R)- 1 . Finally, 

U k u = J-\U k R)(U k R)- 1 . 

This means that the transforms (13.181) induce the following transform of u: 

u i — y U k u. (3.20) 

Now it is easy to check, in view of axial symmetry of A ext , the invariance of L w.r.t. the transforms 
(I3.19p . (13.171) . (I3.20p since p is spherical symmetric. □ 

Recall that e k is the preimage of the basis vector e k w.r.t. the isomorphism ( I2.2ip . The 
Lagrangian L (I3.3P is invariant w.r.t. the spatial rotations (I3.19p . (I3.17p . (I3.20p . In particular, L 
is invariant under the transform groups = e sek . Hence, by the Nother theorem the expression 

M k = M k (X, V) := L v ■ ^p-\s=o (3.21) 

is conserved. 

Definition 3.5 Vector M = (Mi, M2, M3) is called angular momentum of the state (X, V). 
Proposition 3.6 The angular momentum reads 

M = qAq + Iu + J x A E A Bdx + J x A A ext p(x - q)dx. (3.22) 
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Proof For concreteness let us compute with k — 1. Then g s = e sei , and hence 

g s X = (A (e- sil x),e sil A(e- sil x),e sil q,e s£l ). 

Then 

{-he'^x ■ V)A (e- sil x), e^ 1 A(e~ sil x) + e^i-^e'^x ■ V)A{e- sSl x), e^q, ~e x e 



= (eiA Q (x), eiA(x) - (six ■ V)A(x), eiq, ei). 

Note that the last component is the coordinates of ei w.r.t. the basis e~i,e 2 ,e 3 and thus equals 
e\ = (1, 0, 0). Then, since L does not depend on A , 

d 

M 1 = L V - —g s {X)\ s=0 = L A ■ [e x A{x) - (e x x • V)A{x)) + L 4 • {e iq ) + L w • e x 
dx (A ■ (eiA(x) - (e x x • V)A(x)) + VA ■ (eiA(x) - (e x x ■ V)A(x)) 



+q ■ (eiq) + J {eiq) • {A + A ext )p(x - q)dx + Ioj-e 1 + J (e x A (x - q)) ■ [A + A ext }p(x - q)dx 

= {q A g)i + + J {x 2 Af i - x 3 A e 2 xt )p(x - q)dx 

+ J {x 2 A 3 - x 3 A 2 )p{x - q)dx + J (A + VA ) ■ ((0, -A 3 , A 2 ) + (x 3 d 2 - x 2 d 3 )A)dx. (3.23) 

We have to prove that this expression equals to the first component of the RHS of (I3.22p . It 
suffices to prove that the last line equals to the first component of j x A (E A B)dx. Indeed, 
p(x - q) = V • E = V ■ (-Vv4 - A), hence 



(x 2 A 3 - x 3 A 2 )p(x - q)dx = J (x 2 A 3 - x 3 A 2 )(-VA - V 2 A )dx 

= J V{x 2 A 3 -x 3 A 2 ){A + VA )dx. (3.24) 

Then (13.231) transforms to 

J (d 1 {x 2 A 3 - x 3 A 2 )(A 1 + d 1 A ) + x 2 d 2 A 3 {A 2 + d 2 A ) - x 3 d 3 A 2 {A 3 + «9 3 A))) dx 

+ J ({x 3 d 2 - x 2 d 3 )A 1 {A l + d 1 A ) - x 2 d 3 A 2 (A 2 + d 2 A ) + x 3 d 2 A 3 (A 3 + <9 3 A )) dx. (3.25) 
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On the other hand, substitute E = —A — VAq, B = V A A and obtain that the first component 
of J x A (E A B)dx equals 

J x 2 ((d 1 A 3 - 9 3 ^i)(ii + diAo) + (d 2 A 3 - d 3 A 2 )(A 2 + d 2 A ))dx 

- J x 3 ((d 3 A 2 - d 2 A 3 )(A 3 + d 3 A ) + {d x A 2 - a 2 A)(ix + d 1 A ))dx 
which coincides with (I3.25p . The proof is complete. □ 



A Poincare equations 

The derivation of Poincare equations is presented for the convenience of the reader, our exposition 
follows |2]. Poincare has obtained the form of the Hamilton least action principle for Lagrangian 
systems on manifolds [TO] . 

Let V\, . . . , v n be vector fields on a n-dimensional manifold M which are linearly independent 
at every point. Then the commutation relations hold, 

K Vj](g) = Y c Ua)M9), g^M 

where the commutator [i>i,i>j] is defined by 

Kvj}(f) := Vifaif)) - Vj(vi(f)), 

and v(f) is the derivative of a smooth function / on M w.r.t. the vector field v. 

If g(t) is a smooth path in M and / is a smooth function on M, one has g(t) = u>i(t)vi(g(t)) 
and 

j t f(g(t)) = f'(g(t)) ■ g = f(g(t)) • 5>(*)«<GK*)) = 5>(/)w,(f). 

Now consider a variation g(s,t) of the path g(t). Then similarly, 

9 e f{g(e,t)) = Y v j(f) w j( £ ^)^ 
j 

where Wj(e,t) are coordinates of ^(e,t) 6 T g ( E) t)M. Hence 

d e dtf(g(e, t)) = J2J2 v A v iU)) w jUi + Y v i(f) u 'n 

i j i 

d t d £ fe(g(e, t) = Y Y v iMf)) w i u i + u i(/)^> 

3 i i 

where the prime resp. dot stand for the differentiation in e resp. t. However, the differentiations 
in t and e commute, hence we obtain by subtraction 

Y v k(fWk = Y Cij^iWjVkif) + Y v k(f)wk- 

k k ij k 
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Since / is an arbitrary smooth function, we come to the relations 



= ^Cij^iWj + w k . 



(a.i; 



'.I 



Further, let us consider a Lagrangian function L(g,g) on TM. Then L(g,g) can be expressed 
in the variables u: L(g,g) = L(u,g). Let us compute the variation of the corresponding action 
functional taking flA.ip into account: 



— / L(u(e,t),g(e,t))dt = / 



E It"* + V 9 L ■ <i 1 < lf 



E + E 4^) + v ^ • E 



k duJk 



«! k V 



(it dujk 



k 

dL 



dt 



+ E4fc^^T + ^( L ) 



Wk dt. 



The variation should be zero by the Hamilton least action principle, under the boundary value 
conditions 

g(e,ti) = 91, g{e,t 2 ) = g 2 - (A.2) 
Since w k (ti) = w^t^) = by (1A.2j) . we obtain the following Poincare equations: 



d dL 

dt duk 



dL 

"duj 



J2 ( ik UJ ii—+ v k( L )- 



(A.3) 



Remarks 1. If g is expressed in a local map as (<?i, g n ) £ H n , and v k = d 9k , then ( 1A.3I) reduce 
to the standard Euler-Lagrange equations. 

2. If a Lagrangian L does not depend on g, L = L(u) one has 



v k {L) = 0. 



(A.4) 



Indeed, v k (L) = V g L ■ v k (g) = 0. 

3. Suppose M = G is a Lie group, and let i^, fc = 1, n be independent either left-invariant or 
right-invariant vector fields on G. Then c^Ag) are constant: 



(A.5) 



B Commutators of invariant vector fields 

Step 1. By (12.221) the isomorphism (12.211) translates relations (I2.20p to 



[e u e 2 \ = e 3 , [e 3 , e x j = e 2 , [e 2 , e 3 j = ei 
in the sense of matrix commutator. 



(b.i; 
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Step 2. Recall that the right-invariant vector fields v k on SO (3) are defined by right translations 
v k (R) = e k R, where R G 5*0(3) . We should prove (12 .25j) in the sense of the commutators of vector 
fields on the Lie group SO (3). 

Since the fields v k are right-invariant, it suffices to check the relations f !2.25|) at the group unit 
E. Let us compute the derivative of a smooth function / on 50(3) w.r.t. a right-invariant field 
va such that va(E) = A G so(3). In this case v A (R) = AR for R G 5*0(3). Consider a smooth 
path R x {t) G 50(3) such that R^O) = R, i?i(0) = AR. Then 

v A (f)(R) = j t f(Ri(t))\t=o = [f(Ri(t)) • Ri(t)] | t=0 = f'(R) ■ AR. 

In particular, 

V[A,B]{f){E) = f'{E) ■ [A,B], (B.2) 

where [A, B] = AB — BA is the matrix commutator in so(3). 

Now let us compute VA(vB(f))(E) for a right-invariant field vb such that v B (E) = B G so(3), 
v B (R) = BR. Consider a smooth path R 2 {t) G 50(3) such that R 2 (0) = E, R 2 (0) = A. Then 

v A (v B (f))(E) = j t \f{R 2 {t)).BR 2 {t)]\ t=0 = j/(R 2 mt=o-BR 2 m=o+f\R2(tm=o~ t BR 2 (t)\ t=0 



= f'(R 2 (t)) ■ R 2 (t)\ | t =o • BR 2 (t)\ t=0 + [f(R(t)) ■ BR 2 yc )llt=0 

= (f'(E) -A)-B + f{E) ■ BA. 
Then, since the form (f"(E) ■ A) ■ B is symmetric w.r.t. matrices A, B one has 

[v A , v B ](f)(E) = v A (v B (f))(E) - v B (v A (f))(E) = f(E) ■ (BA - AB) = -v [A , B] (f)(E). 

by (IB. 21) . Together with fIB.ip this completes the proof. 
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